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ABSTRACT 



Star formation in the centers of galaxies is thought to yield massive stars with a possibly top-heavy stellar mass distribution. It is 
likely that magnetic fields play a crucial role in the distribution of stellar masses inside star-forming molecular clouds. In this context, 
we explore the effects of magnetic fields, with a typical field strength of 38 yuG, such as in RCW 38, and a field strength of 135 /jG, 
similar to NGC 2024 and the infrared dark cloud G28. 34+0.06, on the initial mass function (IMF) near (< 10 pc) a 10 7 solar mass 
black hole. Using these conditions, we perform a series of numerical simulations with the hydrodynamical code FLASH to elucidate 
the impact of magnetic fields on the IMF and the star-formation efficiency (SFE) emerging from an 800 solar mass cloud. We find that 
the collapse of a gravitationally unstable molecular cloud is slowed down with increasing magnetic field strength and that stars form 
along the field lines. The total number of stars formed during the simulations increases by a factor of 1.5 - 2 with magnetic fields. The 
main component of the IMF has a lognormal shape, with its peak shifted to sub-solar (< 0.3 M G ) masses in the presence of magnetic 
fields, due to a decrease in the accretion rates from the gas reservoir. In addition, we see a top-heavy, nearly flat IMF above ~ 2 solar 
masses, from regions that were supported by magnetic pressure until high masses are reached. We also consider the effects of X-ray 
irradiation if the central black hole is active. X-ray feedback inhibits the formation of sub-solar masses and decreases the SFEs even 
further. Thus, the second contribution is no longer visible. We conclude that magnetic fields potentially change the SFE and the IMF 
both in active and inactive galaxies, and need to be taken into account in such calculations. The presence of a flat component of the 
IMF would be a particularly relevant signature for the importance of magnetic fields, as it is usually not found in simulations. 

Key words. Stars: formation - Stars: mass function - Magnetohydrodynamics (MHD) - X-rays: ISM - Methods: numerical 
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1. Introduction 

Stars are observed to form in molecular cloud fragments 
(clumps) with typical densities of n = 10 4 - 10 6 cm~ 3 . The 
molecular gas in active environments such as galactic cen- 
ters is usually turbulent, FWHM (full-width at half-maximum) 
= ~ 5km/s, and has temperatures ranging from 10 K up to 
a few 1000 K, depending on the local conditions. Molecular 
clouds are observed to have magnetic fields with strengths 
varying from B ~ 10 - few x 1 00 uG (ICrutcher et al.l 119871: 
ICratcherffl999l iBourke et al.ll200lb IWang et al.ll201 ll) . In galac- 
tic centers, it is possible that magnetic field strengths are even 
stronger. From their VLA (Ver y Large Array) observations, 
lYusef-Zadeh & Morrisl dl987albl) estimated using dynamical ar- 
guments that the field strength in our Galactic center could be as 
high as > 1000 uG. A more recent detect i on of mG strengths is 
also reported by iMatsumura et al.l (120121) . ICrocker et all OOlOl) 
show that the minimum magnetic field strength near the Galactic 
center is on the order of ~ 50yuG on 400 pc scales, while they 
infer a minimum field strength of 100yuG using different con- 
straints. Probin g larger scales, rec ent o bservations of the inner 
kpc of galaxies, iBeck et ail (|2005 : ) and iMoss et all d2007l) infer 
magnetic field strengths of B ~ 60 fiG in the central star-forming 
regions of NGC 1097 and NGC 1365, which are amongst the 
strongest fields detected in spiral galaxies. 

Magnetic fields are an important component of the physics 
governing cloud evolution. Yet, it is not fully understood what 



the effects of magnetic fields are during the fragmentation epoch 
of molecular clouds and during star formation. Such early in- 
fluences might also affect the distribution of stellar masses and, 
therefore, the IMF. For a complete understanding of the theory of 
star formation, it is essential to have a complete insight into the 
origin of the IMF. According to most studies over the past two 



important during star and structure formation ( 


Falgarone et al.l 


2001. 2008f Girartetal. 


20091 Wans et alJ 2010 


: Liet al. 2010t 


Schleicher etalJl2010at 


Seifried et al. 201 1; Mouschovias et al. 


201 lb. IWana et al.l d2012l) note, for example, that the shape of 



the filaments of the infrared dark cloud G28.34+0.06 is likely 
affected by magnetic field lines and speculate that the forma- 
tion of the filamentary system has been governed by the in- 
terplay between a strong magnetic fi eld, self-gravity, and tur- 
bulence. Similarly, iGirart et all d2009l) find that the collapse of 
the hot molecular core of the high-mass star-forming region 
G3 1.4 1+0. 31 is controlled by the magnet ic field, resulting in its 
hourglass shape. ISchleicher et al.l d2010al) show that during the 
formation of the first stars, the magnetic field amplification ow- 
ing to the small scale dynamo action will amplify the initial weak 
magnetic field to such an extent that it becomes dynamically rel- 
evant in star formation. Their models show that magnetic fields 
can change the fragmentation properties and also affect the gas 
accretion rates o nto the protostars, bu t only if the resolution is 
sufficiently high dFederrath et alj|20Tlh . Theoretically as well as 
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observationally, magnetic fields are thus shown to play a critical 
role in the evolution and fragmentation of molecular clouds. 

Magnetic fields are generally supposed to slow down col- 
lapse because of the added (Lorentz) force, tension, and pressure 
to the system. A common misconception is that magnetic fields 
should also increase the characteristic mass of stars in a cluster, 
because they tend to increase the effective Jeans ma ss. This 
is not necessarily true. According to the findings of iLi et alj 
(l2OT0h . the magnetic fields rather decrease the characteristic 
stellar mass in a cluster. These authors note that star formation 
occurs along the field lines where dense gas accumulates, 
which increases the density and decreases the Jeans mass, and 
therefore decreases the characteristic mass of the IMF. While 
together with many other magnetic field studies, direct turbulent 
fragmentation is found to be suppressed dHosking & Whitworthl 
12004 iPrice & Bate! 120081: iHennebelle et al.l 1201 ll) . One needs 
to explore further how these results depend on the chosen 
initial conditions and the physics of the simulations. Such 
magnetic field effects could also affect the IMF in star-forming 
regions near massive black holes. Star-forming clouds in these 
environments are, however, subject to strong gravitational 
forces and feedback effects such as irradiation by X-rays 
in active galactic nuclei (AGN), and if inside massive star- 
forming reg i ons, t o increased UV fluxes and cosmic ray rates 
(IWada et al] l2009t iHocuk & Spaansl 120101: Ivan der Werf et al.l 



2010t Meiierink et al.l 1201 U iPapadopoulos et a" 



Hocuk & Spaansl (1201 ll) showed that the impact of radiative 



feedback on the IMF is significant. The IMF that arises from an 
800 M star-forming cloud which is being irradiated by X-rays, 
cosmic rays, and UV photons near a black hole is a top-heavy 
one. 

In this paper, we perform a similar study, where we do not 
investigate stellar feedback, but focus on the effects of magnetic 
fields in the inner 10 pc of galaxies. For this, we simulate model 
clouds with typical magnetic field strengths of B = 38 p.G and 
B = 135yuG, and a model cloud without magnetic fields. By 
comparing these models, we evaluate the impact of magnetic 
fields. 

This paper is organized as follows: In section|2] we introduce 
the numerical code FLASH, construe the implementation of the 
magnetohydrodynamics (MHD) module, and describe the initial 
conditions. In section [3] we give a detailed description on the 
cloud's shape and morphology and how these are affected by 
the black hole's gravity. In section [4] we present our results on 
the effects of magnetic fields for the evolution of the cloud and 
to its star formation. We then look at the resulting IMFs and 
the SFEs from our models and analyze them. We also consider 
radiative feedback in the case that an active black hole irradiates 
the model cloud with X-rays. In section|5] we discuss our results 
and present our final conclusions. 

2. Numerical model and simulation setup 

The simulations in this work have been performed using 
the AMR (adaptive mesh refinement) hydrod ynamical code 
FLASH4 dFrvxell et al.ll2000l: iDubev et al.ll2009h . For modelling 
the magnetic fields, we make use of the recently fr eely available 
'Unsp lit Staggered Mesh' MHD scheme (USM, iLee & Deanel 
120091) . The USM scheme adopts a dimensionally unsplit in- 
tegration on a staggered grid, for the multidimensional MHD 
formulation, based on a finite-volume, higher-order Godunov 
method. This scheme u ses the constrained transport method 
dEvans&Hawlevlll988l) to enforce divergence-free constraints 
of magnetic fields. The USM scheme has shown good per- 



formance in comparison to other current state-of- the-art MHD 
solver s as can be seen from the review paper by iRritsuk et al.l 
(1201 ll) . For solving the Euler equations, we choo se to use the 
HLLD Riemann solver dMivoshi & Kusanol 120051) . which, for 
supersonic MHD turbulence simulations, is among the most ro- 
bust and accurate solvers available, while still staying less diffu- 
sive than other popular solvers like 'Roe' and 'H LLC'. We a lso 
adopt a monotonized central (mc) slope limiter (lLeerlll977b to 
further enhance solution accuracy and stability. 

The hydrodynamical equations are solved for ideal, com- 
pressible MHD conditions. The mass (continuity), momentum, 
energy, and induction equation written in conservative form, 
coupled with the Poisson equation for gravity, are defined as 



dp 

f t + V ■ (pv) = 0, 

^ + V ■ (pvv - BB) + VP = pg, 

at 

^+V-[(E + P)v-(B-v)B]=pg-v, 
at 



» + V.(,B 

ot 



Bv) = 0, 



(1) 



(2) 



(3) 



(4) 



where, p is the fluid mass density, v is the fluid velocity vector, g 
is the gravitational acceleration vector, B is the magnetic field, P 
is the total pressure, and E is the total specific energy. The total 
specific energy, E, and the total pressure, P, are defined as 



r. 1 2 r. B 2 

E = tPV +pE int + — , 
z on 



P = Pth + 
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and the Poisson equation is given by 
V 2 = 47rGp => g = -V0. 



(5) 



(6) 



(7) 



In these equations, Ej nt is the internal energy per unit mass, and 
is the gravitational potential. 

2.1. Particles 

In order to follow individual protostars, high-density gas is con- 
verted into Lagrangian (sink) particles, which are then treated 
as an N-body problem. The sink particles experience forces and 
may themselves contribute to the dynamics of the system. These 
are coupled to the gas through the mesh via gravity. Sink parti- 
cles also have their own motion independent of the grid, so ad- 
ditional equations of motion and forces come into play. Particles 
are advanced using the second order in time accurate, variable 
time step Leapfrog method. The forces acting on the particles 
are handled through 



m i"~jT _ Fiong.i + Fshort.i, 



(8) 



where Fi ongj j represents the sum of all long-range forces (cou- 
pling all particles) acting on the i-th particle, this is handled 
through the mesh, and F s hoit,i represents the sum of all short- 
range forces (coupling only neighboring particles) acting on the 
particle, which is handled directly between the particles. Only 
the gravitational force is implemented for the long and short 
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range interactions. Regular particles are embedded in FLASH. 
Sink particles, on the other hand, are an adapted version of reg- 
ular particles. The difference between them is that sink particles 
can grow in mass by accreting gas and by merging with other 
particles, however, they cannot lose mass or fragment. These 
particles represent compact objects, in our case protostars. 

Sink particles are created when there is locally an irreversible 
collapse. The mass of a sink particle at creation depends on 
the threshold density, which is determined by the Jeans mass. 
This mass also indirectly depends on the maximum grid resolu- 
tion, because the size of the local region is based on a choice. 
Our sink particles typically have a minimum mass of 0.2 so- 
lar masses. This will possibly decrease with higher resolution. 
A sink particle starts to accrete matter immediately after its 
creation. Accretion accounts for the largest portion of the fi- 
nal star mass. This contribution is estimated to be around 85%. 
The mass distribution of stars therefore strongly depends on 
accretion rather than their mass at creation. The sink particle 
algorithm used in this study is further explained in detail in 
iHocuk & Spaansl d2010h. This algori t hm ad opts the same meth- 
ods as described by iRrumholz et all d2004l) and iFederrath et alj 
(I2010h . 

22. Radiation transport and the EOS 

For our primary cloud models (models A, B, and C), we use 
isothermal conditi ons. This is done by using an eq uation of state 
of the form P oc p (IHocuk & Spaansll2010iT201 ll) and fixing the 
gas temperature at 10 K. For our cloud models near an active 
black hole, there is also the presence of X-ray radiation. The 
transport of X-rays and cooling radi ation is handled through a 
radiative transport code (XDR code, Meiierink & Spaansll2005t 
iMeiierink etaUl2007t ISpaans & Meiierinkll2008l) . which is cou- 
pled to FLASH at every time step. For computational speed and 
numerical efficiency, we make a large set of pre-computed tables 
with the XDR code and read it into FLASH during initialization. 
The XDR code incorporates all the heating processes (photo- 
ionization, yielding non-thermal electrons; FUV pumping, fol- 
lowed by collisional de-excitation), cooling processes (atomic 
fine-structure and semi-forbidden lines), and molecular transi- 
tions (CO, H2, H2O, OH, and CH). Cosmic rays, dust-gas cou- 
pling, and secondary photons from X-rays, like internal UV, are 
considered as well. The exa ct implementation of ra diative trans- 
fer is explained in detail in (IHocuk & Spaansll201 ll) . 



2.3. Resolution 

The computational domain encompasses a range of 24 pc in each 
direction, which has outflow boundaries in terms of space and 
isolated boundaries in terms of gravity. Our highest resolution 
is defined by a refinement with a maximum number of 4096 3 
cells, which yields a minimum scale length of 5.7 X 10~ 3 pc, or 
about 10 3 AU. We refine the regions according to a Jeans length 
criterion, where the Jeans length is calculated as 



(9) 



Here, c s is the sound speed, which can be formulated as c^ = 
yP/p or cj; = yT for an ideal gas. y is a parameter which de- 
pends on the equation of state (EOS). For a polytropic EOS, the 
pressure scales as P = Kp y . Ay- 1 in this case represents 
an isothermal state. Whenever the Jeans length of any cell in a 



block of 16 3 cells would fall below 20 times the cell size, re- 
finement is initiated for the whole block. When A\ is resolved by 
over 50 times the cell size, de-refinement is initiated in order to 
save computational time and memory. 

2.4. Initial conditions 



Similar to the performed studies in lHocuk & Spaansl d2010t) and 
IHocuk & Spaansl (120111) . we simulate a gravitationally unstable 
spherical molecular cloud at dbh = 10 pc from a 10 7 M black 
hole and follow its collapse. The initial cloud has a number den- 
sity of 10 5 cirT 3 with a mean molecular mass p = 2.3 and a 
radius of 0.33 pc. These conditions create a cloud of 800 so- 
lar masses. The rest of the 24 3 pc sized simulation box has a 
number density of 10 cirT 3 giving a total mass of 8000 M . The 
cloud is in a Keplerian orbit around the black hole with an or- 
bital velocity of VGM/dbh = 66.7 km/s. This introduces some 
gravitational shear into the cloud. The shearing velocity across 
the cloud is given by 



AVshear 



(10) 



where G is the gravitational constant, r c j ouc j is the cloud radius, 
and M is the total enclosed mass within dbh, which is dominated 
by the black hole, i.e., M ^ Mbh- In this equation, dbh » fcioud- 
The imposed shear across the model cloud is on the order of 
2.2 km/s. This affects the cloud morphologically and adds tur- 
bulence to the system. The shearing turbulence is comparable 
to, but slightly lower than, the initial turbulence of 5 km/s. The 
shearing time, t s h ea r = 2r c i ou d/Av S hear, is about three times larger 
than the initial free-fall time of the cloud, where tff is given by 



t ff 



3tt 
32Gp 



10 s yr. 



(ID 



Thus, a gravitationally bound molecular cloud of 10 5 cirT 3 is 
able to exist in such an environment. As mentioned before, 
we initiate the molecular cloud with turbulent conditions that 
are typical for an active environment, FWHM = 5 km/s, or 
0"turb =2.1 km/s (iPerez-Beaupuits et al.1 120091) . and apply this 
over all scales with a power spectrum of P(k) o c k~ 4 , follow- 
ing the empirical laws for compressible fluid s (lLarsoniri981l : 
iMvers & Gammielll999t iHever & Bruntll2004l) This scaling is 
also k nown as Burgers turbulence dBurgersI 1939t lB ec & Khaninl 
120071) . Besides the injected energy of gravitational shear, we do 
not drive the turbulence in this work. 

For the initial magnetic field, we associate our model cloud 
with the observed s tar-forming region RCW 38 surround ed by 
a molecular cloud dBourke et all 1200 ll IWolk et al.1 120081) with 
a magnetic field strength of B = 38 pG. We adopt this value 
and initialize a uniform magnetic field strength throughout the 
domain. The orientation of the initial magnetic field vectors is 
chosen to be a uniform field with equal strength along each 
spatial direction, i.e., the field lines make an angle of 45 de- 
grees with each axis. To study and compare against molecular 
clouds with higher magnetic field strengths, we compute another 
model with B = 135 pG. The higher field strength is compara- 
ble to molecular clouds such as NGC 2024 and the infrared dark 
cloud G28.34+0.06. These molecular clouds have similar den- 
sity (10 5 cm 4 ) and size (fewxO.l pc) as the model cloud. 
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Table 1. Model parameters 



Model 


B-field strength 


Magnetic flux O 


M/O 


Magnetic pressure 


Tad 


X-ray flux 




ImG] 


[Gauss cnT 2 ] 


[wrt critical] 


[dyne cnT 2 ] 


[yr] 


[erg s" 1 cm- 2 ] 


Model A 



















Model B 


38 


1.20 X 10 32 


21.76 


5.75x10-" 


5.20 x 10 7 





Model C 


135 


4.24 x 10 32 


6.12 


7.25x10-'° 


6.57 x 10 s 





Model D 


38 


1.20 x 10 32 


21.76 


5.75x10-" 


5.20 x 10 7 


160 


Model E 


135 


4.24 x 10 32 


6.12 


7.25x10-'° 


6.57 x 10 s 


160 



Note: The primary cloud models are models A, B, and C. 



We assume that the magnetic flux is frozen to the gas. Flux- 
freezing holds even for primordial clouds as long as the mag- 
netic force is smaller than the gravitational force during col- 
lapse (Maki & Susa 2004, 2007), and should become even better 
with increasing ionization degree due to cosmic rays or radiative 
backgrounds. In order to be able to collapse, the ratio of gaseous 
mass over the magnet i c flux needs to be larger than a critical 
value. IFalgarone et al.l d2008i) derived the mean mass-to-flux ra- 
tio A for a sample of 14 clumps of massive star formation using 
their measured line of sight magnetic field strengths. They found 
this ratio to be A 6 ± 0.5, where A is defined as 



A=-/{-\ 



(12) 



Here, <D is the m agnetic flux, (M/<I>) cr j t is the c ritical mass-to- 
flux ratio, and c$ dMouschovias & Spitzei 19761) is a factor that 
is exp ected to be about 0.12 — (27 T T 1 dMouschovias & Spitzed 
fl976t iNakano & Nakamural \l9W. iTomisakaet al.l 119881) . We 
have adopted the value of (27r)-' for c$. For the model clouds 
in this study, given our fixed clump mass of 800 M and field 
strengths of 38 fiG and 135 //G, we have mass-to-flux ratios 
of A = 21.8 and 6.1, respectively. This means that both model 
clouds are supercritical (A > 1, i.e., prone to gravitational col- 
lapse). The model parameters are listed in TableQ] 

Considering that the magnetic pressure is P mag = B 2 /87r, the 
initial pressure exerted by the Lorentz force inside the cloud be- 
comes P mag = 5.75 x 10- 11 dyne cirr 2 for model B (38 fiG) 
and 7.25 x 10 10 dyne cm 2 for model C (135 fiG). Since the 
thermal pressure for an isothermal temperature of 10 K is P t h = 
1.38 x 10~ 10 dyne crrr 2 , these pressures yield a f3 = P t h/P m ag 
2.40 for model B and /3 = 0.19 for model C. 

The temperature of the primary cloud models remain isother- 
mal at 10 K. The isothermal sound speed of the cloud is c s = 0.19 
km s _1 in these runs. We also run the simulations B and C for 
the case of an active black hole and name them models D and 
E. An active black hole will produce strong X-rays (1 - 100 
KeV) in its accretion disk. A 10 7 M black hole accreting at 
10% Eddington radiates with a bolometric flux of ~10 4 erg s 1 
cirr 2 at 10 pc distance. Most (90%) of the energy is emitted 
at optical and UV wavelengths, but these wavelengths will be at- 
tenuated along the line of sight by large columns of gas and dust. 
Such environments are typical of obscured AGN, as for med in 
(U)LIRGS, with obscuring columns of 1Q 22 -10 23 - 5 cm 2 dAaltol 
120051: iPerez-Beaupuits et al.ll2Q07t [Loenen et al1l2008l) . X-rays, 
however, can penetrate large columns of gas and dust and dom- 
inate the t hermal balance up to d istances of 300 pc from the 
black hole (ISchleicher et al.l2010bl) . The expected (1-100 KeV) 
X-ray flux at 10 pc from the black hole is ~160 erg s cm 2 
dHocuk & Spaansl 120101) . The cloud models with X-rays ob- 
tain their tempera tures from the XDR co de, which can go up 
to ~ 5 x 10 3 K dHocuk & Soaansl 120101) . Fig. Q] displays the 



temperature-density phase diagram obtained from one of the X- 
ray models. Under such conditions, the sound speed will rise to 



1.E-03 



5.E-01 



1.E+01 



* 3 



o 2 




log n (cm 3 ) 

Fig. 1. Temperature-density phase diagram of an X-ray irradi- 
ated simulation (model D) at t = 2 tff. The color in this figure 
represents the amount of mass (in M ) lying in a contoured re- 
gion. 



a maximum of 5 km/s in the hottest parts of the X-ray irradi- 
ated cloud. A temperature of a fewxlO 3 K results in a thermal 
pressure of P t h - 8 x 10~ 8 dyne cirr 2 , which is 2 - 3 orders of 
magnitude higher than the magnetic pressure, /? 1 10 - 1400. 

Ambipolar diffusion according to the conditions of our 
model is not expected to play an important role for our calcu- 
lations, because it operates on a much longer timescale than the 
free-fall time of tff — 10 5 yr. The ambipolar diffusion timescale 



r ad = 3 x 10 6 yr 



\10 4 cm- 3 / {': 



B 



-2 



fcloud 



30yuG/ \0.1pc 



(13) 



is on the order of r a d = 5.2 - 65.7 x 10 7 yr for the two mag- 
netized cloud models with magnetic field strengths of 38 and 
135 yuG. The ambipolar diffusion can work in clouds with lower 
ionization degree or at smaller scales like compact cores inside 
the cloud. In such a small area, r c i ou d is much smaller and the 
magnetic fiel d is much stronge r (but t he density will be larger). 
According to iNakamura & Lil d2005l) . the ambipolar diffusion 
can be accelerated by the turbulent compression and therefore 
may play a role in the evolution of cores out of which stars form, 
in particular in turbulent clouds. For our X-ray models, the ion- 
ization degree will be higher which will make ambipolar diffu- 
sion less efficient. 
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Model-A 



Model-B 



Model-C 





1.0 - 




0.8 - 


o 
Q. 




CD 


0.6 - 


N 




tfl 


0.4 - 




0.2 - 








■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I 


■ 1 1 1 1 1 1 1 1 , 




dX . 

dZ 

(dX+dY)/2 








■ ... i ■ ... i ■ ... i ■ ... i ■ ... i ... ■ 



0.0 0.5 



1.0 1.5 2.0 
time (t„=10 5 yr) 



2.5 3.0 




1.0 1.5 2.0 
time (t„=10 5 yr) 



1.0 1.5 2.0 
time (t„=10 5 yr) 



Fig. 2. Spatial sizes of the cloud. The size of the cloud in each spatial direction is plotted against time in tff . The sizes in the direction 
of rotation are denoted by dX and dY, and the height is given by dZ. The dashed line is a representation of the average disk size. 



3. The collapsing cloud 

3.1. Cloud morphology 

A molecular cloud that is gravitationally unstable will collapse 
in about one free-fall time, which is about 10 5 years in all our 
simulations. By virtue of geometrical effects for an object in or- 
bit around a gravitating body and due to the tidal shear imposed 
by the gravity of the black hole, any symmetrical collapse is dis- 
rupted. Clouds in a strong gravity field will collapse in one di- 
rection first and quickly form a disk. Furthermore, the existing 
gravitational shear generates turbulence and pressure in the two 
spatial directions of the orbital plane, while the third spatial di- 
rection (Z-direction) experiences no additional pressure and can 
proceed with a near free-fall collapse. In Fig. [2] we plot the size 
of the cloud from the X, Y, and Z perspective against time. From 
the cloud sizes in this figure, we can see that a disk forms in the 
plane of the orbit within one free-fall time. The increase in den- 
sity due to this asymmetric collapse creates even more pressure 
along the orbital plane of the cloud. In the same figure we also 
see that this pressure causes a slight expansion of the cloud in 
the XY-plane. 

The cloud itself is selected by taking the gas densities above 
10 5 ctrT 3 and the sizes are selected for densities above the mean 
cloud density. The collapse in the Z-direction continues unim- 
peded for each model, but slows down between 1 - 1 .5 tff. The 
disk thickness stabilizes around 0.1 pc in all runs. We can see 
that the shape of the cloud in the X and Y directions stays 
roughly the same. The expansion of the disk, which can be in- 
ferred from the dashed line ((dX + dY)/2), is explained by the 
gravitational shear stretching the cloud along the orbital plane 
(the XY-plane) in the direction of rotation and by the collapse in 
the Z-direction which increases the pressure in the same plane. 
We see that disk formation is not affected by the presence of 
magnetic fields of the chosen magnitudes in this work. 

The magnetic field lines align parallel to the disk in less than 
a free-fall time as the disk forms. The alignment of the mag- 
netic field with the disk is because the gravity is more important 
than the magnetic force. The vector components along the XY- 
plane get enhanced by the increase in density in the course of 
gravitational collapse. Moreover, because of tidal shear from the 
black hole, magnetic field components parallel to the disk will 
form rather quickly and only the parallel components are am- 
plified when the field is compressed by the flow. The field lines 
are completely aligned with the disk around 1 — 1.5 tff, while, at 
the same time, the field lines in the orbital plane shape like an 
hourglass. An image to this effect is shown in Fig. [3] 




Fig. 3. Magnetic field vectors versus disk orientation. The mag- 
netic field lines align with the disk between ~ 1.0-1.5 tff. The 
image shows a volume plot of the density from an XZ perspec- 
tive (top) and an XY perspective (bottom) at t = 1 .5tff . The mag- 
netic field vectors are displayed as white arrows, where the field 
strength B typically ranges from 1-1 000 fiG. The X-axis is given 
at the bottom in units of parsec. 



3.2. Turbulence 

The model cloud develops turbulence during its evolution as a 
result of the strong gravity field of the black hole. To demon- 
strate this, we show radial profiles of the rms velocity and the 
rms density of the cloud for three different time frames in Fig. 
2] All models start with a cloud that is turbulent with a Mach 
number of M — 9 and a number density of 10 5 cirT 3 . From 
the changing velocity field in the different snapshots we can see 
that the cloud gets more turbulent in time with Mach numbers 
that grow to around M = 50. This means that the turbulence 
increases by a factor of 5 - 6 within the time frame of 3 tff. 

We also notice that the rms velocities are initially lower at 
small radii and higher at large radii. This can be inferred from the 
rms velocities displayed in the top left panel of Fig. [4] Since we 
take the mean of the velocity at each radius to compute the rms, 
it is appropriate to state that the radius is analogous to the scale 
length. At t = 2 tff, the turbulence at smaller scales is increased 
by the decay of turbulence from larger scales. 

We find that the rms velocities are lower with increas- 
ing magnetic field strength. This is especially true at smaller 
radii. Given that the density is the highest close to the center 
of the cloud, the magnetic field strength is also greater there. 
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Fig. 4. Radial profiles of the cloud. The rms velocities and the rms densities of models A, B, and C are plotted against cloud radius. 
The radius originates in the center of the cloud. The rms velocities are obtained in the frame of reference of the cloud orbiting the 
black hole. From left to right, the time frames of t = 1,2, and 3 tff are displayed. 



Pseudocolor 

Vat: pat iicleSolat mass 

0.24 0.29 0.35 0.42 




Max: 0.51 
Min: 0.24 



Vector 
Var: Magvec 
0.51 2.3e-06 8.0e-05 0.00016 0.00024 0.00031 




Max: 0.00031 
Min: 2.3e-06 





Fig. 5. Star formation along the field lines. Sink particles (col- 
ored spheres) are aligned with the magnetic field (white arrows) 
for model B at t = 1.1 tff, i.e., just after the first sinks have 
formed. The masses of sink particles are given in solar mass units 
and the magnetic field strength is given in Gauss. The axes are 
given in parsec. 



4. Results 

4.1. Star formation along magnetic field lines 

We see that most of the sink particles form along the field lines 
and, especially, where the field lines are the strongest, which is 
generally the central part of the cloud. The first sink particles 
are created at t ^ 1 tff around densities of n > 10 6 cm -3 , with a 
thermal Jeans mass of Mj < 0.5 M Q . The magnetic field strength 
has increased by a factor of 10 - 12 from the initial value at the 
onset of star formation. Fig.[5]illustrates the sink particle forma- 
tion along the magnetic field lines just after the first sink parti- 
cles have for med. These results are in accordance with earlier 
findings (e.g., Nakamura & Lill2"008h iLi et al.ll2010h . The reason 
why star formation is occurring along the field lines is explained 
by the directionality of the Lorentz force. The magnetic field in- 
creases the plasma pressure, by an amount B 2 /87r, in directions 
perpendicular to the magnetic field, and decreases the plasma 
pressure, by the same amount, in the parallel direction. Thus, 
the magnetic field gives rise to a magnetic pressure, B 2 /87r, act- 
ing perpendicular to field lines, and a magnetic tension, B 2 /87r, 
acting along field lines. This allows the gas density to build up 
along the field lines until it becomes gravitationally unstable. 



4.2. Star-formation efficiencies 



Additionally, the rms densities are also smaller with magnetic 
fields. These lower rms values indicate that there is less cloud 
fragmentation with increasing magnetic field strength. In all 
cases, the mean density is also lower at the same moment in time 
with increasing magnetic field strength. From the reduced den- 
sity build-up, we can understand that magnetic fields success- 
fully slow down collapse. However, the difference between the 
three models in the rms density and the mean density decreases 
after t = 1.5 t ff . 



Sink particle formation: We display the SFEs of the primary 
models in Fig. [6] The two lines in each panel of this figure sketch 
the time evolution of the total sink particle mass (solid line) and 
the total amount of accreted mass onto sink particles (dashed 
line) with respect to the total gas mass. Examining the SFE of 
model A, we see that it fo llows a typical behavior that is also 
seen in previous studies (IClark et al.l 120081 : iHocuk & Spaansl 
l20lTh . The onset of star formation for model A is at t = 0.85 tff 
and the SFE at t = 3 t ff is 10.7%. The final SFEs for models B 
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Fig. 6. SFEs of the primary models. The ratio of the total sink particle mass over the total mass is plotted against time in tff as the 
solid line. The red dashed line portrays the total accreted mass in relationship to the total mass. In the frames from left to right, the 
initial magnetic field strength increases from 0, 38, to 135 //G. The total number of sink particles formed during the run is given in 
the upper right corner of each panel. 



and C are, on the other hand, somewhat lower. Both magnetic 
field runs enjoy a final SFE of ~8%. Typical sink particle masses 
at the time of their creation is around 0.3 - 0.8 M Q for model 
A, 0.2 - 0.5 M for model B, and 0.2 - 0.3 M Q for model C. In 
model A, we obtain a total number of sink particles, denoted as 
np, of 215 as measured at t = 3 tff. In the presence of magnetic 
fields, the number of sink particles increases drastically by a 
factor of 1.5 - 2, e.g. np = 435 for model B. We presume this 
increase to reflect the changes induced in the global dynamics, 
in particular with respect to cloud morphology, gas reservoir, 
and turbulence. When going to an even larger field strength 
(model C), the number of sink particles reduces to 330, due to 
the decreased SFE and the effect of magnetic pressure. 

Sink particle accretion: The relatively uniform initial condi- 
tions in the cloud prompt it to form the bulk of its sink parti- 
cles in about a free-fall time. The sink particles that have al- 
ready formed disrupt the formation of new sink particles in their 
vicinity by accreting the material around them and by disturbing 
the ge neral process of sink particle formation dHocuk & Spaansl 
201 1). Because of this, star formation is virtually terminated be- 
tween 1.3 — 1.7 tff. 

We note that the SFE curves flatten out just after the onset of 
star formation around 1.3 tff and pick up again above 2 tff. The 
delay that the magnetic fields cause to star formation clearly has 
a strong effect on the star-formation history. Magnetic fields re- 
duce the gas flow rates into the center of the cloud by magnetic 
pressure and also lower the initial stellar masses at their creation. 
The lower initial masses as well as the reduced gas supply both 
lead to reduced accretion rates in models B and C. We also see 
this in the SFEs, which is always smaller in the presence of mag- 
netic fields. As accretion scales with the square of the particle 
mass, having twice more particles with twice less mass means 
that the total accretion onto protostars will be reduced by a fac- 
tor of 2, thus reducing the typical mass scale of the stars. We can 
see from Fig. [6] that the amount of accreted mass is smaller for 
magnetic field runs in the first 10 5 years after the onset of star 
formation. On the other hand, this also means that the remain- 
ing gas reservoir will maintain a high mass for a longer period. 
Around t = 2 tff, we see a jump in the SFEs in the magnetic 
field runs, which may suggest that there are regions in the cloud 
that collapse at different epochs. Rather than a second burst of 
star formation, we find that this is mainly due to accretion. Sink 
particles start to accrete more mass in the magnetic field runs be- 
cause of the higher remaining density (see Fig. but also due 
to mergers happening in the center of the cloud, where most of 



the stars lie. Mergers become more likely when the masses of 
stars are higher. Since the newly merged stars from this epoch 
are more massive, they tend to accrete more gas as well. 

4.3. The initial mass functions 

The IMFs are constructed by counting the number of stars within 
a mass range of dM = 0.133 in log units. For this, we chose a 
mass range of 0.2 M Q to 20 M and divided it into 16 bins. In or- 
der to make the selection independent of epoch, we count all the 
sink particles between one and three free-fall times, with a time 
resolution of l/60th tff, and take the average. Alternatively, we 
have also analyzed the IMFs by selecting the time frames where 
all the models have converted the same amount of gas mass into 
sink particles. Either way, we get the same results. The IMFs 
of the primary models are shown in Fig. [7] I n these images , the 
green dashed lines indicate the Salpeter IMF (ISalpeterl 19551) and 
the bl ue dot-dashed lines represent the Chabrier IMF dChabrierl 
120031) . The red solid and the purple triple-dotted lines are the 
best power-law and the best lognormal fits to the data. The IMFs 
derived from the simulations with (models B and C) or without 
(model A) magnetic fields differ significantly. The run without 
magnetic fields resembles the conditions of the Milky Way with 
an IMF that nicely follows a lognormal shape and has a power- 
law slope of T = -1.38 above 1 M . This is also in agreement 
with ea rlier numerical studies of the IMF for Milky Way con- 
dition s dEimegreenll 19971; iBonnell et al]|2006t iHocuk & Spaansl 
l2010h . The characteristic mass is slightly larger than expected, 
i.e., m c hai ~ 0.9 M , where the commonly found value in the 
Milky Way is around 0.3 - 0.5 M . This could be the result of 
the strong gravity field in our models. 

For the runs with magnetic fields, the IMF is mainly pop- 
ulated by sub-solar mass stars < 0.3 M . While power-law fits 
below ~ 1 M yield slopes of T = -2.75 and -3.27, the IMF 
seems no longer well represented by a single power-law or a 
lognormal function. Although there is still a log-normal com- 
ponent, its peak is shifted to smaller mass scales, due to the 
decreased accretion rates discussed in the previous subsection. 
For model C, the sink particle distribution indeed seems to re- 
semble an exponential form, as the low-mass tail of the IMF is 
no longer resolved. Besides the lognormal component, magnetic 
fields tend to create an additional flat component of the IMF 
above 1-2 solar masses, as a result of the additional pressure. 
This second component is a local effect from the magnetic fields, 
which stabilizes against local collapse until higher densities are 
reached. The effect seems to occur only locally where the mag- 
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Fig. 7. The IMFs of the primary models. The images display the time-averaged IMFs between 1 and 3 free-fall times, where 
tff = 10 5 yr. From left to right, the magnetic field strength increases from 0, 38 to 135 jt/G. In each image, for comparison purposes, 
the Salpeter IMF (green dashed) and the Chabrier IMF (blue dot-dashed) are displayed (with arbitrary values along the Y-axis). Two 
best fits are applied to the data, a power-law fit and a lognormal fit, and are shown as purple and red lines. With the exception of the 
lognormal fit, the slopes above the characteristic mass are given in the upper left corner. 



netic field strength is particularly high, while most of the formed 
sinks contribute to the log-normal component in regions domi- 
nated by turbulence. This behavior makes sense, considering that 
the mean mass to flux ratio is still larger than one in all simula- 
tions. These findings markedly coincide with the bottom-heavy 
IMFs observed in massive cluster e l liptica l galaxies at high red- 
shift (Ivan Dokkum & Conrovl[2010ll201 ll) . 

4.4. An active black hole 

We tested two additional models in which the black hole was 
considered to be active. The accreting black hole produces 
strong X-rays that impinge on the model cloud with a flux of 
160 erg s _1 cm 2 . The two X-ray models, D and E, have the 
same parameters as models B and C (Table Q]) aside from be- 
ing irradiated by X-rays. The temperatures in the X-ray models 
range from 5500 K at the irradiated face of the cloud to 10 K 
in the most shielded regions. We show the sizes, an IMF, and an 
SEE in Figs. [8] [9] and[TU] 

If we look at the cloud morphology and compare models D 
and E (Fig. [8]) against B and C (Fig. O, we see that the sizes are 
more compact in each direction. This is a direct consequence of 
X-ray irradiation. X-ray heating increases the external pressure 
by raising the temperature of the low-density ISM and particu- 
larly the cloud surface. This aids collapse. At the same time it 
also decreases SFEs (Fig. ITOb. because (1) X-ray heating evap- 
orates the outer layers of the cloud, which leads to mass loss, 
and (2) X-rays increase the temperature inside the cloud in re- 
gions where the radiation can penetrate, and by doing so, the 
Jeans mass is raised. We find that the SFE of model D is only 
1% of the total mass whereas model E shows no star forma- 
tion at all. For model E, the combined pressure of the magnetic 
field and the increased internal thermal pressure due to X-rays, 
where P mag +th = 2.71 x 10~ 8 dynecirT 2 at t = 1 tff, is halting 
collapse completely, with P grav = 2.33 x 10~ 8 dynecirT 2 in the 
cloud core, and a slowly evaporating, quiescent disk forms. In 
both X-ray models, the cloud loses a significant portion of its 
mass to the ISM, over 60% at t = 1 tff, especially at densities of 
n < 10 5 cm 4 

If we examine the IMF of model D in Fig. [9] we see that the 
IMF has a power-law slope F of -1.30 and a characteristic mass 
of 0.37 M Q . The characteristic mass in this model peaks at the 
same mass as in model B. The global effect of magnetic fields 
on accretion has the same implications for model D. However, 



the number of stars is considerably reduced (by a factor 10). We 
see that X-ray irradiation has strongly quenched star formation 
in this cloud. The existence of a temperature gradient throughout 
the cloud is the underlying reason. The formation of low-mass 
stars is hindered on the irradiated side of the cloud and can only 
exist in the shielded, cold regions. The more massive stars are af- 
fected to a lesser degree. Stars with M > 2 M have not formed 
due to the largely decreased SFE and a flat high-mass contribu- 
tion is not seen. These effects altogether change the IMF back to 
a log normal shape w ith a Salpeter slope, in a manner conforming 
with iKroupal (120021) . 

If typical clouds have ~ 800 M , as considered here, it im- 
plies that predominantly the log-normal component of the IMF 
will appear, as the overall efficiency is too low for the formation 
of high-mass stars. In larger clouds, one may however still ex- 
pect to see the flat component of the IMF as well, provided that 
the gas reservoir is sufficiently high. 

5. Summary and Conclusions 

We presented a series of numerical simulations to study the im- 
pact of magnetic fields on the IMF in star-forming clouds near a 
black hole. To this end, we created a magnetically supercritical, 
800 M cloud at 10 pc from a 10 7 M o black hole and followed 
its evolution. One model was simulated without magnetic fields 
as a fixed basis for comparison (model A). This model had sim- 
ilar conditions as the galactic center of the Milky Way. We find 
that model A results in a universal IMF with typical Salpeter 
slope of T = -1.38. The characteristic mass ofm c h ar ~O.9M is 
somewhat on the large side. We attribute this to the presence of 
a strong gravitational field from the black hole. To test the im- 
pact of magnetic fields, we initiated the cloud with two different 
magnetic field strengths, 38 and 135 fiG (models B and C). All 
three models were performed under isothermal conditions. 

In all the models, the cloud contracts into a disk within a 
free-fall time, but the collapse is slower with magnetic fields. 
We find that the magnetic field lines align parallel to the disk 
as the disk forms and arrange to an hourglass shape. The align- 
ment of the magnetic field with the disk is because gravity is 
more important than the magnetic force. The vector components 
along the XY-plane get enhanced by the increase in density in 
the course of gravitational collapse, during which the turbulence 
aids the deformation of the field lines across the molecular cloud. 
We see less turbulence with increasing magnetic field strength 
and find that direct turbulent fragmentation is reduced. The first 
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Fig. 8. Spatial sizes of the cloud. The size of the cloud in each spatial direction is plotted against time in tff . The sizes in the direction 
of rotation are denoted by dX and dY, and the height is given by dZ. The dashed line is a representation of the average disk size. 
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Fig. 9. The IMF of the X-ray model D. The image displays 
the time-averaged IMF between 1 and 3 free-fall times, where 
tff = 10 5 yr. In this image, the Salpeter IMF (green dashed) 
and the Chabrier IMF (blue dot-dashed) are displayed (with 
arbitrary values along the Y-axis). Two best fits are applied to 
the data, a power-law fit and a lognormal fit, and are shown as 
purple and red lines. The slopes above the characteristic mass 
are given in the upper left corner. 



Fig. 10. The SFE of the X-ray model D. The ratios of the to- 
tal sink particle mass (solid line) and the total accreted mass 
(dashed line) with respect to the total mass is plotted against 
time in tff. The total number of sink particles formed during 
the run is given in the upper right corner. Note that the plotted 
range along the Y-axis differs from the range in Fig. [6] 



sink particles form along the field lines, but dynamic interac- 
tions disperse them quickly. The number of sink particles that 
form in the magnetic field runs are considerably higher (factor 
of 2 for B = 38 yuG) with smaller masses. The reason is that 
the slower transport of gas into the center is reducing individual 
stellar masses and their accretion rates and thus more sink par- 
ticles can form, while the overall change in cloud morphology, 
gas reservoir, and turbulence seems to favor enhanced fragmen- 
tation. We also notice that there is a jump in the SFEs around 2 
tff in the magnetic field models (Fig. |6). Sink particles start to 
accrete more mass at this epoch, because the density increases 
and mergers are happening in the center of the cloud. 

The IMFs are severely affected by magnetic fields. We find 
that the IMFs are very bottom heavy with increasing magnetic 
field strength, but have an additional flat contribution at the 
high-mass end. Our best fits show very steep power-law slopes 
below ~ 1 M H , consistent with h igh redshift observations 
(Ivan Dokkum & Conrovf2010ll201 ll) . We also find a decreasing 
trend in the characteristic mass in the presence of magnetic 
fields, as a result of the decreasing accretion rates. The respec- 
tive characteristic masses for models B and C are 0.35 M Q and 
0.20 M H . A d ecrease in characteristic mass is also found by 
iLi et al.l (120101) in their numerical simulations. So overall, both 
the decreasing stellar masses and the additional flat component 
in the IMF are in the end due to the magnetic pressure. The 
first is more because of its global effects, which decreases the 



accretion rate, and the second is a local effect, which stabilizes 
against local collapse until a higher mass is reached. 

We summarize our main conclusions for the isothermal magnetic 
field simulations as follows. When there are magnetic fields: 

• Collapse is slowed down. 

• Star formation occurs along the field lines. 

• Fragmentation is reduced. 

• The number of stars increases by a factor of 1.5 - 2. 

• The SFE decreases by ~ 20%. 

• The characteristic mass shifts to smaller mass scales. 

• The lognormal IMF is very bottom heavy (T < -2.75). 

• The IMF has an additional flat high-mass component. 

It is e xpected that magnetic fields allow the formation of 
filam ents dFalgarone et al.l 120011: iBanerjee et ail 120091: ILi et al.l 
120101) . We do not see any clear evidence of filaments in our mag- 
netic field simulations. Our different initial conditions might be 
affecting their formation. Because of the strong gravity of the 
nearby black hole, the cloud quickly collapses into a disk. The 
injected turbulence by the tidal shear imposed by the black hole 
is also influencing the shape of the cloud and likely disallowing 
filament formation. 

The initially chosen magnetic field orientation could be 
of importance for the results. The orientation of the initial 
magnetic field vectors in this work is a uniform field with 
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equal strength along each axis. This means that the field lines 
make an angle of 45 degrees with each axis. We have tested 
the importance of the initial orientation by running one case 
with a different orientation, i.e., a uniform magnetic field where 
the field lines are perpendicular to the orbital plane. There 
are some differences in results, such as the number of formed 
sink particles, but we find that the main conclusions (cloud 
morphology, IMFs, and SFEs) still stand. An in-depth study on 
the importance of the initial magnetic field orientation is beyond 
the scope of this work. 

X-ray feedback by an AGN: We have also simulated a realistic 
scenario in which the black hole is active and accreting at 10% 
Eddington. A large fraction of the AGN's primary output is 
obscured by interstellar gas and dust close to the accretion disc, 
but the most energetic wavebands, like X-rays, can penetrate 
large columns of gas and dust and irradiate the cloud. We have 
used an X-ray flux of 160 erg s _1 cnT 2 , which is typical for 
obscured AGN environments. This setup is performed with the 
two aforementioned initial magnetic field strengths of 38 and 
135 nG. 

We summarize our main conclusions for the X-ray models as 
follows. When there are X-rays: 

• Collapse is aided by the increased external pressure. 

• The cloud slowly evaporates and the Jeans mass is increased. 

• The number of stars is reduced by an order of magnitude. 

• The overall SFE is strongly inhibited. 

• The IMF has a lognormal shape with T = — 1 .30. 

• For higher-mass clouds, one may expect an additional flat 
component of the IMF, as seen in the isothermal runs. 

Our results tell us that in the presence of magnetic fields, 
star formation in galactic centers experiences a mode where low- 
mass stars are preferred, but that X-ray feedback reduces the 
SFE. For 800 M Q clouds as considered here, X-rays make the 
power-law slope of the IMF flatter and allow a typical Chabrier 
IMF to emerge. We expect significant changes in the IMF in 
galactic centers with inactive massive black holes both due to the 
shifted lognormal as well as the additional flat component, but a 
decreased SFE and a quenched low-mass star formation in AGN. 
As observations may have limitations in detecting the low-mass 
stars, a flat (top-heavy) IMF can be expected in the presence of 
magnetic fields since, especially at the low-mass end of the IMF, 
completeness i s often an issue, but significant progresses are be- 
ing made (e.g. jDa Rio et al 1 12009b . Evidence for such a compo- 
nent would thus provide strong evidence for the importance of 
magnetic fields during star formation, as hydro dynamical mod- 
els ge nerally yield lognormal distributions (e.g. lNakamura & LJ 
120081) . Increasingly accurate measurements of the IMF in the 
center of our own galaxy may thus provide a relevant pathway to 
probe their implications (|Dib et al. 2007; Espinoza et al.ll2009t 
iHuBmann et alj2~012tlLiermann et aLl2012l) . 
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